Nanocrystal Formation in Si Implanted Thin Si0 2 Layers under the Influence of an Absorbing 

Interface 
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Kinetic 3D lattice Monte Carlo studies are presented on Si nanocrystal (NC) formation by phase separation in 
1 keV Si + implanted thin Si0 2 films. The simulation start from Si depth profiles calculated using the dynamic, 
high-fluence binary collision code TRIDYN. From the initial Si supersaturation, NCs are found to form either by 
nucleation, growth and Ostwald ripening at low Si concentrations. Or at higher concentrations, non-spherical, 
elongated Si structures form by spinodal decomposition, which spheroidize by interface minimization during 
longer annealing. In both cases, the close Si0 2 /Si interface is a strong sink for diffusing Si atoms. The NCs 
align above a thin NC free oxide layer at the Si0 2 /Si interface. Hence, the width of this zone denuded of 
NCs has just the right thickness for NC charging by direct electron tunneling, which is crucial for non- volatile 
memory applications. Moreover, the competition of Ostwald ripening and Si loss to the interface leads at low Si 
concentrations (nucleation regime) to a constant width of the denuded zone and a constant mean NC size over a 
long period of annealing. 
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I. INTRODUCTION 

Nanocrystals (NCs) embedded in a matrix are known to 
form by phase separation of a supersaturated solid solution. 
At low impurity concentrations, nucleation (jj and conserva- 
tive Ostwald ripening [2] will lead to an coarsening ensemble 
of precipitates during annealing . According Lifshitz, Sly- 
ozov, and Wagner OllJh me mean particle size increases ac- 
cording a power law in time. At higher concentrations, phase 
separation takes place by spinodal decomposition |5J. A van- 
ishing nucleation barrier |6] implies that small initial con- 
centration fluctuations amplify during the annealing. Non- 
spherical structures precipitate, which form a connected net- 
work of precipitates above the percolation threshold 0]. In 
the presence of an absorbing (reflecting) interface phase sepa- 
ration is considerably altered. For spinodal decomposition, it 
is known that the interface directs the wave vector of concen- 
tration fluctuations |8]. Hence, a layered structure of precipi- 
tates forms close to the interface. Recent studies also revealed 
an influence on nucleation, which can either be enhanced or 
depressed by the absorbing/reflecting interface |9]. 

Phase separation in thin films is of particular interest for 
the fabrication of non-volatile memories based on MOS tran- 
sistors containing semiconductor NCs in their gate oxide 11(11 . 
Using especially 1 keV Si + implantation into 8 nm thick Si0 2 
on (001)Si, NCs were formed a few nanometers above the 
Si/Si0 2 interface lITlll . which is in direct tunneling distance. 
Thus, these NCs can be charged by direct electron tunneling 
from the Si, which is a prerequisite for a high endurance and 
low operation voltages of NC-memories [12]. Nevertheless, 
the influence of the absorbing Si0 2 /Si interface on Si NC for- 
mation is not well understood and requires further investiga- 
tions. For that aim, kinetic 3D lattice Monte Carlo (KMC) 
simulations of Si phase separation from Si0 2 were reported 
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Figure 1: TRIDYN depth profiles for 1 keV Si + implantation into 
Si0 2 (denoted by open and bold symbols). The profiles were fitted 
to Gaussian distributions (Table|Q shown as solid curves. 

recentlv fl3ll and shall be discussed in more detail in this pa- 
per. As input data, Si depth profiles for 1 keV Si + implanta- 
tion into Si0 2 were taken from dynamic, high-fluence binary 
collision simulations conducted with TRIDYN 1 1 
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II. BINARY COLLISION SIMULATION OF HIGH 
FLUENCE SI+ IMPLANTATION 

The TRIDYN depth profiles of implanted Si into of Si0 2 
are shown in Fig.^for 1 keV energy. Similar to the widely 
applied binary collision simulation code TRIM 1 15], TRIDYN 
is used to calculate ion range profiles in amorphous targets 
but includes additionally dynamic target changes as induced 
by sputtering, swelling and ion beam mixing. The input pa- 
rameters being required include the displacement and surface 
binding energies of the target atoms. For both, Si and O, the 
displacement energy was assumed to be 8 eV. The surface 
binding energies of Si and O are assumed to vary linearly with 
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the surface composition that they balance the enthalpies of 
sublimation and decomposition of Si and Si0 2 , respectively 

GEO. 

The Si depth profile broadens with increasing fluence due 
to swelling, sputtering and mixing (Fig.^ Additionally, pref- 
erential sputtering has lead to a Si enrichment at the target 
surface. Nevertheless, the profile shape remains Gaussian in 
the fluence range considered here. The parameters of Gaus- 
sian distributions fitted to the TRIDYN profiles of Fig.^ are 
summarized in Table|I] which also includes the calculated tar- 
get swelling. 



III. THE KINETIC MONTE CARLO METHOD 

The implanted Si separates from the Si0 2 matrix during 
subsequent annealing, i.e. precipitates are formed. In general, 
this process of phase separation consists of a huge number 
of elementary events (like bond breaking, diffusional jumps 
of atoms, chemical reactions etc.) that occur in random se- 
quence. Their simulation must take into account both ener- 
getic and statistical aspects, which is here suitably conducted 
by kinetic 3D lattice Monte Carlo (KMC) simulations. The 
KMC method applied is discussed elsewhere in detail 1 18], a 
short overview will be given instead. 

The kinetics of Si atoms is described in a solid host ma- 



Table I: Parameters of Gaussian distributions fitted to TRIDYN pro- 
files of Fig.Q Additionally, the target swelling rf swe u ing due to incor- 
porated Si ions is given. 
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trix (Si0 2 ), which is the background or "system's vacuum". 
Thereby, an underlying fee lattice has been assumed, which 
together with hep is the most isotropic lattice. Within this 
host dissolved Si diffuses and can form precipitates. (Crystal- 
lites obey the same four-fold symmetry than in the diamond 
lattice. The lattice spacing were chosen such that the correct 
atomic Si density is obtained.) Applying the classical lattice 
gas model with attractive Si-Si interaction, the energetics is 
determined by the nearest-neighbor (NN) Ising model. The 
kinetics of the system is governed by the Kawasaki particle 
exchange dynamics using the Metropolis algorithm 11911 . Sta- 
tistically, each Si atom is allowed to jump from an initial site 
i to an empty neighboring site / during one Monte Carlo step 
(MCS), which is the time unit of the KMC simulation. The 
transition probability P.* for one atom to jump is given by 

_ JV e ~ £ A "/>«,■, 

P 'f = 1^-1 -\ E A+( n i-"f) E B\P 

[T e L A • >' B v , tif < n t , 

where z^ 1 denotes the attempt frequency, n i * represents the 

number of NN bonds at the two sites and j3 = (kT)~ l has its 
usual meaning. E A is the diffusion barrier for Si in Si0 2 and 
E B the Si-Si bond strength. In principle, the bond strength E B 
can be determined from the solubility of Si in Si0 2 via the 
detailed balance of Si attachment/detachment at the Si/Si0 2 
interface. (It has to be taken into account that the coordi- 
nation number in the fee lattice of the simulation is 12 in- 
stead of 4 in the case of the Si lattice.) In kinetic Monte 
Carlo simulations it is convenient to renormalize the transition 
probability of the most probable event to one. For E A — > 
each diffusional jump is allowed, which defines the temper- 
ature dependent time scale of a Monte Carlo step (MCS), 
T = T Q e\p{E A j5}. Then, the dimensionless transition proba- 
bility reads as P^ = min(l, exp | — (n i — n^)E B fi j). Energeti- 
cally unfavored events (« ( > iij) are allowed, nevertheless, to 
occur by thermal activation according the Boltzmann statis- 
tics. 
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Figure 3: Snapshots of KMC simulations of phase separation in 8 nm thick Si0 2 on (001) Si implanted with Si. The initial depth profiles were 
taken according Fig.Qfor fluences of 5 x 10 15 cm~ 2 (a) and 1 x 10 16 cm~ 2 (b). Phase separation proceeds via spinodal decomposition (a,b) 
and is accompanied by percolation at the highest fluence (b). The 15 nm scale applies only in the lower right corner. 



In principle, the scale of time and temperature of the sim- 
ulation is given by the experimental parameters of the system 
under consideration. However, the diffusivity and the solu- 
bility of Si in Si0 2 are largely unknown, thus reduced units 
are applied for the time - the Monte Carlo step - and for the 
temperature - the reduced temperature £gj3. Nevertheless, 
the path of systems's evolution towards equilibrium predicted 
by KMC simulations may improve the process understanding 
substantially. 



IV. KMC SIMULATIONS OF PHASE SEPARATION 

In Fig. |2 and |3 snapshots of the KMC simulations of the 
phase separation are shown for E B fi = 2. Thereby, the TRI- 
DYN profiles of Fig. [2 were used as initial Si distribution. 
The size of simulation volume is 56 x 56 x Az nm 3 , where Az 
is 8 nm plus swelling due to implantation according Tab. U 
For the sake of clarity, only a quarter of the simulation cell 
is shown in Fig. [2] and [3] besides of the lower right corner 
in Fig. |51 Periodic boundary conditions were applied in the 
plane, while at the surface a reflecting boundary conditions 
was applied. The simulation cell (Si0 2 ) borders on fixed (001) 
layers of Si atoms accounting for Si attachment/detachment at 
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Figure 4: Evolution of mean NC shape measured by 8 in the nu- 
cleation regime (2 x 10 15 cm~ 2 ) and for spinodal decomposition 
(5 x 10 15 cm -2 ). The highest fluence (1 x 10 16 cm~ 2 ) is omitted due 
to percolation. 



the Si/Si0 2 interface, which is an effective sink for diffusing 
Si atoms. 
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Two regimes of phase separation are found by KMC sim- 
ulations. A "nucleation and growth" regime is observed 
for 2 x 10 15 Si + cm~ 2 . Si NCs nucleate and grow further 
at the expense of Si supersaturation (Fig. [2] 100 kMCS). 
Later on, NCs grow by Ostwald ripening and finally dissolve 
due to Si loss to the SiO z /Si interface, respectively. Above 
5 x 10 15 Si + cm~ 2 a "spinodal decomposition" regime is iden- 
tified, where due to the vanishing nucleation barrier | fj non- 
spherical, elongated Si structures form (Fig.[3]a, 100 kMCS). 
At even higher Si concentrations (1 x 10 16 Si + cm~ 2 ), above 
the percolation threshold, these structures become laterally 
connected and extend over several tens of nanometers (Fig. [3] 
b, 100 kMCS). 

As seen in Fig. [2] and [3] nucleation and spinodal decom- 
position lead to different initial shapes of precipitates, which 
can be accessed by defining a shape parameter 8 1 18] (analo- 
gously to the moment of inertia), 



It measures the mean squared distance of atoms from the cen- 
ter of mass of the NC relative to that of a sphere with ra- 
dius R, which assu mes that all N atoms within a NC form 
a sphere with R = y/ (3N)/(4n). Thus, the shape parameter 
8 should tend to one for a spherical NC. Non-spherical struc- 
tures should obey 8 > 1. The time evolution of the shape 
parameter 8 is shown in Fig. [4] In the nucleation regime, the 
NC shape deviates little from the spherical form (8 « 1.2) at 
the beginning, while it approaches quickly a constant value 
of 5 « 1.1 during further annealing. (Small NC may have a 
8 slightly larger than one due to their facetted shape.) For 
spinodal decomposition 8 the formed structures are consider- 
ably non-spherical (8 > 1.3) initially. Below the percolation 
threshold, non-spherical Si precipitates evolve into spherical 
NCs (5^1) due to interface minimization, which can hardly 
be distinguished from NCs formed by nucleation (see Fig. [2] 
and|3(a), 1000 kMCS). 

The close Si0 2 /Si interface acts as an effective sink for Si 
diffusing within the Si0 2 in both regimes, which results in 
a zone denuded of NCs formed at the interface. However, a 
more detailed consideration reveals differences between nu- 
cleation and growth and spinodal decomposition. As shown 
in Fig. [5] for the nucleation regime (2 x 10 15 cm~ 2 ) the width 
of the denuded zone is constant over a long period of anneal- 
ing. A similar behavior is observed for the mean NC size 
(Fig- ID- The competition of NC dissolution due to Si loss to 
the interface and Ostwald ripening leads to a constant mean 
NC diameter. 

On the other hand, the mean precipitate size is larger in 
the spinodal decomposition regime (Fig. [3}. According the 
Gibbs-Thomson relation, the equilibrium Si solubility in Si0 2 
at the precipitate's interface is lower than in the nucleation 
regime. Thus, the diffusional flux of Si form the precipi- 
tates towards the Si/Si0 2 interface is greatly reduced. Conse- 
quently, the NC diameter is still increasing. Within the anneal- 
ing period shown in Fig. [6] Ostwald ripening is more effective 
then the Si loss to the Si0 2 interface. Moreover, interface min- 
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Figure 5: Width of the zone denuded of NCs between the Si0 2 /Si 
interface and the NCs. 
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Figure 6: Evolution of the mean NC diameter and the NC density 
during annealing for fluences of 2 x 10 15 cm -2 and 5 x 10 15 cm~ 2 . 
The highest fluence (1 x 10 16 cm~ 2 ) is omitted due to percolation. 



imization of non-spherical precipitates leads to a narrowing of 
the denuded zone (Fig. [3}- 

In both regimes, NCs form at high density as shown in 
Fig- El if the percolation threshold is not reached. In the ini- 
tial stage, the NC density is slightly higher in the nucleation 
regime, while during longer annealing the NC dissolve faster 
than in the spinodal regime. 

For a fluence of 2 x 10 15 Si + cm~ 2 (nucleation regime), the 
evolution of mean NC size and density is compared for phase 
separation in 8 nm and 20 nm thick Si0 2 (Fig. 0. For both 
case, the same ion energy of 1 keV and, thus, the same ini- 
tial Si depth profiles were used. Up to an annealing time of 
10 6 MCS, the NC size and density obeys the same power law 
for both oxide thickness (seen as straight line in the Log-Log 
plot of Fig. 0. The Si/Si0 2 interface does not influence the 
early NC evolution. Later on (> 1000 kMCS), NCs start to 
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Figure 7: Comparison of mean NC size and density for phase sepa- 
ration in 8 nm and 20 nm thick Si0 2 layers implanted at a Si fluence 

of 2 15 cm" 2 . 



dissolve in the 8 nm thick Si0 2 layer. Thus, the mean NC size 
becomes constant and NCs dissolve faster than by Ostwald 
ripening, which is seen as kink in the plot of the NC density 
after 2 x 10 6 MCS. The evolution of the NC size and density 
is not affected by the far Si0 2 /Si interface in the case of the 
thick (20 nm) Si0 2 . 



Si from Si0 2 . Two major regimes of phase separation, nu- 
cleation and growth, and spinodal decomposition were ob- 
served. At high concentrations, spinodal decomposition is ac- 
companied by percolation. The regimes of phase separation 
determines the initial shape of the precipitates. Nucleation 
and growth leads to almost spherical NCs, while spinodal de- 
composition results in non-spherical Si structures, which may 
spheroidize during long lasting annealing. 

The performed process simulations predict that for non- 
volatile memory applications Si NCs should be preferably 
synthesized in the nucleation regime. Then, NCs form at 
high density (> 10 12 cm -2 ). The influence of the absorbing 
Si/Si0 2 interface leads to NC dissolution and a zone denuded 
of NC forms at the interface. This denuded zone has just the 
right thickness of 2 . . .4 nm in order to act a barrier for NC 
charging by direct electron tunneling from the Si. Moreover in 
the nucleation regime, the width of this barrier is found to be 
constant over a long period of annealing. A similar behavior 
holds for the mean NC size. For spinodal decomposition, the 
situation is different. The width of the denuded zone and the 
mean NC size varies during annealing (Fig. [5] and [6}. Above 
the percolation threshold, the formed elongated structures are 
connected laterally, and, therefore, would behave like a float- 
ing gate in a conventional MOS transistor. Charge brought to 
this network of precipitates would spread and one oxide de- 
fect could lead to a complete discharging of the floating gate, 
which could not occur in the case of laterally well-isolated Si 
NCs. 



V. SUMMARY 
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